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Abstract Force due to the self-gravity of the disc in 
the vertical direction is considered to study its possi- 
ble effects on the structure of a magnetized advection- 
dominated accretion disc. We present steady-sate self 
similar solutions for the dynamical structure of such a 
type of the accretion flows. Our solutions imply re- 
duced thickness of the disc because of the self-gravity. 
It also imply that the thickness of the disc will increase 
by adding the magnetic field strength. 
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1 INTRODUCTION 

There has been rapidly progress over the past three 
decades towards a better understanding of the accretion 
processes in astrophysics, in particular accretion discs 
around compact objects or even black holes (see re- 
views by Narayan, Mahadevan & Quataert 1998; Kato, 
Fukue & Mineshige 2008). Black hole accretion is one 
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of the most important ingredient when considering as- 
trophysical scenarios of galaxy /quasar formation. Cur- 
rent belief is that there are two sites where an accretion 
disc around a black hole can be found: In close binary 
systems called X-ray binaries (XBs), and at the cen- 
ter of the galaxies. The birth of modern accretion disc 
theory is traditionally attributed to the original model 
presented by Shakura & Sanyev (1973). This standard 
geometrically thin, optically thick accretion disc model 
(SSD) can successfully explain most of observational 
features in active galactic nuclei (AGNs) and X-ray bi- 
naries. In the standard thin disc model, the motion 
of matter in the accretion disc is nearly Kcplerian, and 
the viscous heat in the disc is radiated away locally. An 
alternative accretion disc model, namely, the advection- 
dominated accretion flows (ADAFs), was suggested for 
the black holes accreting at very low rates (Ichimaru 
1977, Narayan & Yi 1994). 

In some of the scenarios of structure formation in as- 
trophysics, in particular those related to the formation 
of the stars or galaxies, self-gravity of the system may 
play a vital role. Formation of an accretion disc, as one 
stage of the structure formation, is an important part of 
any theory in this field. However, in the standard accre- 
tion disc model, the effect of self-gravity in the vertical 
or radial direction of the disc is neglected for simplicity 
and the disc is supported in the vertical direction only 
by the thermal pressure. Although in some of the ac- 
creting systems it is a reasonable assumption, there are 
situations, in which one can hardly neglect self-gravity 
of the disc itself. According to the current theories, 
self-gravity of the disc not only can contribute in gen- 
erating turbulence inside the disc but its force in the 
radial and the vertical directions would modify the an- 
gular velocity and and the vertical scale- hight, respec- 
tively Duschl et al. 2000. Therefore, it is not an easy 
task to include self-gravity of a disc in a self-consistent 
approach. In fact, the theory of self-gravitating ac- 
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cretion discs is less developed. Early numerical work 
of self-gravitating accretion discs began with N-body 
modeling (Cassen & Moosman 1981; Toomley; Cassen 
& Stein-Camcron 1991). Although great progress has 
been made in recent years in increasingly sophisticated 
numerical accretion disc simulations, simple analytic 
disc models still are the only accessible way of mak- 
ing direct link between the theory and observations. 
The theoretical treatment can estimate the spectra and 
other observational features of the accretion powered 
objects. To solve the nonlinear equations of the self- 
gravitating accretion discs, sometimes the technique of 
self-similarity is useful. Some of the astrophysical sys- 
tems often attain self-similar limits for a wide range of 
the initial conditions. Moreover, the self-similar prop- 
erties allows us to investigate properties of the solu- 
tions in arbitrary details, without any of the associated 
difficulties of numerical magneto-hydrodynamics. Sev- 
eral classes of the self-similar solutions for the structure 
of self-gravitating accretion discs have been studied by 
now (e.g., Paczynski 1978; Mincshigc & Umemura 1996; 
Shadmchri 2004) . Some of them are presented for the 
standard accretion discs (e.g., Paczynski 1978), while 
the others are describing ADAFs including self-gravity 
(e.g., Mincshigc & Umemura 1996; Shadmchri 2004). 
Abramowicz et al. 1984 examined a number of issues re- 
lated disc's self-gravity in a steady state accretion disc, 
especially in regard to the effect of disc self-gravity on 
the disc topology and disc dimensions. Bu, Yuan & 
Xie (2009) have proposed a self-similar solution of a 
magnetized ADAFs without self-gravity. Mineshige & 
Umemura (1996) presented a set of self-similar solu- 
tions for steady-state structure of ADAFs considering 
gravitational force of the disc in the radial direction. 
However, such solutions are hardly applicable because 
the gravitational force of the disc due to its self-gravity 
in ADAfs is generally smaller than the radial force of 
the central object. However, the gravitational force of 
the disc in the vertical direction can be a fraction of the 
vertical component of gravitational force of the central 
object, at least at the outer part of ADAFs. Then, as 
we will show, the structure of the disc would modify be- 
cause of the self-gravity of the disc itself. We think, our 
solutions are applicable to the outer parts of ADAFs in 
X-ray binaries. 

Although we did not apply our solution to a specific 
astronomical system, we believe the solutions are appli- 
cable to the discs with an inner hot part surrounded by 
an outer cold part. In the most cases, one can hardly 
recognize the transition region. For example, the outer 
regions of the X-ray binaries are SSD, though their in- 
ner parts resemble to the hot accretion flows (e.g., Esin 
et al. 1997). Our solutions are suitable for a region 



where the disc is self-gravitating but the advection of 
the dissipated energy is not negligible. Also, the discs 
of the low- luminosity AGNs (e.g., Ho 2008) may have a 
similar situation and they do not show clear transition 
from SSD to ADAF. The solutions are also applicable 
to Slim discs, and thus have its potential applications 
in ULXs and NLSls (Mineshige et al. 2000; Watarai 
et al. 2000). The next step for future studies is to 
relax our simplifying assumptions, in particular simi- 
larity method, and solve the relevant equations for one 
of the mentioned systems and in doing so, our similarity 
solutions will guide us. 

2 The Basic Equations 

We are interested in analyzing the dynamical behav- 
ior of a magnetized ADAF when the force due to the 
self-gravity of the disc in the vertical direction is con- 
sidered. But self-gravitational force in the radial di- 
rection is neglected because ADAFs do not extend ra- 
dially so much. We consider only toroidal component 
of the magnetic field. We suppose that the gaseous 
disc is rotating around a compact object of mass M*. 
Thus, for a steady axisymmctric accretion flow, i.e. 
d/dt = d/d(j) = 0, we can write the standard equations 
in the cylindrical coordinates (r, <f>, z). We vertically 
integrate the flow equations and so, all the physical 
variables become only functions of the radial distance 
r. We also neglect the relativistic effects and Newto- 
nian gravity in radial direction is considered. The disc is 
supposed to be turbulent and possesses an effective tur- 
bulent viscosity v. As for the energy conservation, we 
assume the generated energy due to viscosity dissipa- 
tion is balanced by the radiation and advection cooling 
(e.g., Narayan & Yi 1994). 

The equation of continuity gives, 

-^-{rY,v r ) = 2pH (1) 
r dr 

where v r is the accretion velocity (i.e., v r < 0) and p 
denotes the mass-lose rate per unit volume due to the 
wind or outflow, H is the disc half-thickness, and E is 
the surface density. We can also write E = 2pH. 
The equation of motion in the radial direction is 

^dVr = < _ GJW, _ 1 d. (Ec 2j _ C \ _ J_±^ c 2 ) 

r dr r r 2 E dr s r 2E dr A ' 

(2) 

where v v , c s and ca are the rotational velocity of the 
disc, sound and Alfvcn velocities of the gas, respec- 
tively. Sound speed and the Alfven velocity are defined 
as c 2 = Pga S /p and c\ = B 2 /4irp = 2p mag / p where 
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p gas and Pmag are the gas and the magnetic pressures, 
respectively. 

The vertically integrated angular momentum equa- 
tion becomes 



rEv r — (rv v ) = — (r 3 i/£— ) 
ar ar ar 



(3) 



where Q = v^/r is the angular velocity. Also v repre- 
sents the kinematic viscosity coefficient and we assume 
(Shakura & Sunyaev f973), 

v = ac s H (4) 
where a is a nondimensional parameter less than unity. 

The goal of this study is to investigate the effect 
of self-gravity of the disc on its structure. One can 
estimate the importance of self-gravity of the disc by 
comparing the contributions to the local gravitational 
acceleration in the vertical direction by both the central 
object and the disc itself. The vertical acceleration at 
the disc surface due to its self-gravity is 2irGT, and by 
the central object is GAI^h/r 3 . Thus, force due to the 
self-gravity of the disc in vertical direction is dominated 
if: 



M d 7rr 2 S 1 h 
M7 ~ Ml > 2~7 



(3) 



where Md is the mass enclose the disc within a radius 
r. For ADAFs the typical values of h/r is around 1. 
Since the enclosed mass Md is an increasing function of 
the radial distance r, the effect of self-gravity becomes 
significant when the disc is thick. Here we consider ke- 
plcrian sclfgravitating (KSG) disks in which selfgravity 
is significant only in the vertical direction and which 
satisfy the constraint (1/2) (ft/r)M» < M d (r) < M*. 
Compared to the standard ADAF solutions (Narayan 
& Yi 1994; Akizuki & Fukue 2006), the KSG disks re- 
quire modification of the equation of hydrostatic sup- 
port in the direction perpendicular to the disk. Thus, 
while in the standard ADAF solutions the local vertical 
pressure gradient is balanced by the z component of the 
gravitational force due to the central object, in the SG 
case we have balance between two local forces, namely 
the pressure force and the gravitational force due to the 
disks local mass. Also in the KSG case, in the radial 
direction centrifugal forces are still balanced by gravity 
from a central mass (Keplerian approximation) . 

For a self-gravitating disc, the hydrostatic equilib- 
rium in the vertical direction yields (e.g., Paczyncki 
1978, Duschl ct al. 2000) 



where P c is the pressure in the central plane (z = 0) 
and G is the gravitational constant. Since details of 
the thermodynamics in the vertical direction are not 
considered in our model, we shall assume the disc to be 
isothermal in the vertical direction. 

For a magnetized Keplerian self-gravitating disc, we 
have 
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where p = ^ is a vertically averaged mass density and 
p _ P^ag _ i / ca \ 2 wnere ^his parameter shows the 
important of magnetic field pressure in comparison to 
the gas pressure. We will take j3 as an input parameter 
of our model and its effect on the structure of the disc 
is studied for different values of /3. Having Equations [6] 
and we can write 



H 



2ttGE 

Also, the energy equation becomes 

EtV dc 2 s + f£_^_/ x _ j^/d^-yi 
7 — 1 dr r dr r dr 



(8) 



(9) 



where 7 is the ratio of specific heats. The advection 
parameter / measures the degree to which the flow is 
advected-dominatcd (Narayan & Yi 1994), and it is sup- 
posed to be constant. 

Induction equation with the field scape can be writ- 
ten as (Akizuki & Fukue 2006) 



— (V r B v ) — B v , 



(10) 



where B v is the field scaping/creating rate due to the 
magnetic instability or the dynamo effect. We can 
rewrite this equation as 



dc A 
dr 



C4 



dv r 
dr 



c 2 A v r 



2cl 



B4 
B d 



2 2 P H 



(11) 



Now, we have a set of ordinary differential equations 
that describe the dynamical behavior of a magnetized 
ADAF flow including self-gravity of the disc. Solu- 
tions of these equations give us the dynamical behavior 
of the disc, which depends on the viscosity, magnetic 
field strength, self-gravity and advection rate of energy 
transport. 



3 Self-Similar Solutions 



P c = ttGS 2 



(6) 



Global behavior of the disc can not be described us- 
ing similarity solutions, because boundary conditions 
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Fig. 1 These plots show behavior of the coefficients 
Ci(— and Ci{^) versus self-gravitating parameter D 
for different values of the advection parameter /. Here, the 
input parameters are p = 0.1, a = 0.1 and 7=1. 



Fig. 2 These plots show behavior of the coefficients 
Ci(— and C2(^ L ) versus self-gravitating parameter D 
for different values of the magnetic field parameter /3. Here, 
the input parameters are / = 0.7, a = 0.1 and 7 = 1. 



are not considered in this method. But, as long as we 
are not interested in the solutions near the boundaries, 
similarity solutions describe correctly, true and useful 
asymptotically behavior of the flow in the intermediate 
regions. 

We assume that the physical quantities can be ex- 
pressed as a power law of the radial distance, i.e. r v , 
where v is determined by substituting the similarity 
solutions into the main equations and solving the re- 
sulting algebraic equations. Therefore, we can write 
similarity solutions as 



v r (r) = —CiOVk(r), 



v v {r) = C 2 v k (r), 



ct(r)=C 3 v 2 k (r), 



(12) 



(13) 



4(r) = 
where 

Vk(r) = 



B 2 v (r) 
4irp(r) 



2/3C 3 v 2 k (r), 



GAI 



(15) 



(16) 



and the parameters C\,Ci and C3 will be determined 
from the main equations. 



In addition, the surface density E is assumed to be 
a form of 



(17) 



where So and s are constants. Then, in order for the 
self similar treatment to be valid, the mass-loss rate 
per unit volume and the escaping rate must have the 
following form, 



P = Par 



-5/2 



B,! 



(18) 



(19) 



where po and Bq are constants. 

Considering hydrostatic equation, we obtain the disc 
half-thickness H as 



(14) H _ C 3 (l + /3) _ C 3 (l 



9 Md 



D 



(20) 



2|| i ) as self-gravitating 



where Md is the mass enclosed in the disc within a 
radius r and is given approximately by Md = irr 2 T,. 
Also, we define parameter D( 
coefficient. 

By substituting the above similarity solutions into 
the continuity and the motion equations and also the 
energy and the induction equations, we obtain the fol- 
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lowing system of dimensionless algebraic equations, 



1. daXoVGMtD 



(21) 



-a 2 Cl + Cl - 1 - [s - 1 + 0(s + 1)] C 3 = 0, (22) 



ct = 0.1 , p = 0.1 , 7=1.0 



f - 0.3 
f - 0.5 
- f = 0.7 



3(8 + !)(! + 0)^/2 



D 



3 > 



H (1 + /3)C 3 
r D 

2 _ 2D , 3-7 1/2 

^-9/(TT^) ( 7^T )Cl 3 ' 



B = 



(3-s) 



aCi(GM Ht ) 1 



~0 + W' 



(23) 
(24) 
(25) 

(26) 



As is easily seen from the above equations, for s = —1/2 
there is no mass loss, while there exist mass loss for 
s > —1/2. On the other hands, the escape and cre- 
ation of magnetic fields are balanced each other for 
s = 3. Although outflow is one of the most important 
progresses in accretion theory,(see Narayan & Yi 1995; 
Blandford & Bcgclman 1999; Stone & Pringle & Begcl- 
man 1999 and also some recent work like Xie & Yuan 
2008; Ohsuga & Mineshige 2011), but in out model we 
chouse self-similar solution in a way that p = and 
B v oc r~ n / 4 ,(s = —1/2). So we ignored the effect of 
wind and outflow on the structure of the disks. 

After some algebraic manipulations, we can reduce 
the above equations to a sixth order algebraic equation 
for Ci, 

A 3 Cf + 3A 2 ECf + (3AE 2 + B 3 )C 2 + E 3 = 0, (27) 

where the coefficients depend on the input parameter 
as 



E 



1 



2D 



( W a (- '-) 

35/3/4 + ^ S-r 



(28) 

(29) 
(30) 
(31) 
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Fig. 3 These plots show behavior of the h/r versus self- 
gravitating parameter D for different values of the advection 
parameter / upper panel and magnetic field parameter, f3, 
lower panel. Here, the input parameters are / = 0.7, a = 
0.1, P = 0.1 and 7 = 1. 



Having C\ from this algebraic equation, the other vari- 
ables (i.e. C*2 and C3) can be determine easily. 

Now, we can do a parameter study considering our 
input parameters. Figure 1 shows how the coefficients 
Ci and C2 depend on the self-gravitation parameter D 
for different values of advection parameter /. Radial 
velocity is determined by Ci which is shown in the up- 
per panel. In ADAFs the radial velocity is generally 
less than free fall velocity on a point mass, but it be- 
comes larger if the advection parameter / is increased 
and such a behavior is consistent with the previous an- 
alytical solutions (e.g., Akizuki & Fukue 2006; Abbassi 
et al. 2008, 2010). Obviously, larger parameter D im- 
plies a more self-gravitating disc and top plot of Figure 
1 shows that the accretion velocity decreases as param- 
eter D increases. Compare to a non-sclf-gravitating 
disc, low D, the velocity of the flow will decrease some 
order of magnctudc. Moreover, as more dissipated en- 
ergy is advected within the flow, we see that the radial 
velocity is more sensitive to the variations of the self- 
gravitating parameter within the range D ■< 5. But 
when the disc becomes strongly self-gravitating (i.e., 
large D), the accretion velocity is more or less indepen- 
dent of the advected energy. Bottom plot of Figure 1 
shows coefficient C2 (i.e., rotational velocity) versus pa- 
rameter D. Here, the rotational velocity is not affected 
by the variation of the parameter D, but it decreases 
when the advection parameter / increases. 
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Fig. 4 These plots show behavior of the coefficients 
Ci(— and C2{-^) versus advection parameter / for dif- 
ferent values of the self-gravitating parameter D. Here, the 
input parameters are /? = 0.1, a — 0.1 and 7 = 1. 

In Figure 2, we assume the fraction of the advected 
energy is / = 0.7 and then behavior of the coefficients 
C\ and C2 versus parameter D are shown for differ- 
ent values of magnetic field parameter /?. For a given 
D, both the radial and rotational velocities increase as 
the disc becomes more magnetized (i.e., larger /?) A 
magnetized disc must rotate faster than a case without 
magnetic field because of the effect of magnetic tension. 
The radial velocity is more sensitive to the variations 
of D when the disc is more magnetized. The rotational 
behavior is the same for different values of D but it 
slightly shifts up when parameter j3 increases. 

In Figure 3, we show vertical thickness of the disc in 
terms of the self-gravitating parameter D for different 
values of / and f3. The disc becomes thinner as the 
parameter D increases. If we neglect force due to the 
self-gravity of the disc in the vertical direction, the disc 
becomes thicker as more energy advected with the flow. 
But our analysis shows for a strongly advected disc even 
not very large values of D lead to a reduction of the disc 
thickness by a factor of two. For a given D, by increas- 
ing the parameters / and (3 the vertical thickness in- 
creases which means that advection and magnetic field 
will cause the disc becomes thicker. 

Figure 4 shows how the coefficients C\ and C2 de- 
pend on the advection parameter / for different val- 
ues of self- gravitating parameter D. Radial velocity is 
determined by C\ which is shown in the upper panel. 
In the ADAFs solutions usually by adding advection 



Fig. 5 These plots show behavior of the coefficients 
Ci(— and Ci(^) versus advection parameter / for dif- 
ferent values of the magnetic field parameter /3. Here, the 
input parameters are D = 0.2, a — 0.1 and 7 = 1. 

parameter radial velocity will increases while the rota- 
tional velocity decrease ( Akizuki & Fukuc 2006). In- 
fluence of the self-gravitation of the disc will lead to 
decrease of radial velocity while the rotational velocity, 
C2, is not affected. In Figures 5 and 6 the influence of 
magnetic field strength, /3, on the behavior of CiS and 
vertical thickness of the flow were plotted respectively. 
By adding (3, which indicates the role of magnetic field 
in the dynamic of accretion discs, we will see that the 
radial flow increase as well as rotational velocities and 
vertical thickness increase. On the other hand radial 
and toroidal velocity increase when the toroidal mag- 
netic field becomes large. This is due to the magnetic 
tension term, which dominates the magnetic pressure 
term in the radial momentum equation that assist the 
radial in-fall motion. 



4 Conclusions 

In this paper, we studied an accretion disc in the ad- 
vection dominated regime considering purely toroidal 
magnetic field and the force due to the self-gravity of 
the disc in the vertical direction. A set of similarity so- 
lutions was presented for such a configuration. Our so- 
lutions reduce to the previous analytical solutions (e.g., 
Akizuki & Fukue 2006; Fukue 2004; and Abbassi et al 
2008) when self-gravity is neglected. Some approxima- 
tions were made in order to simplify the main equations. 
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Fig. 6 These plots show behavior of the h/r versus self- 
gravitating parameter D for different values of the self- 
gravitating parameter D upper panel and magnetic field 
parameter, /3, lower panel. Here, the input parameters are 
P = 0.1, a = 0.1, D = 0.2 and 7 = 1. 



we believe the solutions are are presented, are appli- 
cable to the discs with an inner hot part surrounded 
by an outer cold part where the self-gravity plays an 
important role and and advection of the energy is not 
negligible. The discs of the low-luminosity AGNs may 
have a similar situation and they do not show clear 
transition from SSD to ADAF. The solutions are also 
applicable to Slim discs, and thus have its potential 
applications in ULXs and NLSls. 

It would be interesting if we add the effect of the 
wind and outflow in our future investigation to find 
out how it would change solutions. The next step for 
future studies is to relax our simplifying assumptions, 
in particular similarity method, and solve the relevant 
equations for one of the mentioned systems and in doing 
so, our similarity solutions will guide us. 
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We assume an axially symmetric and static disc with 
a-prescription for the viscosity. We also ignored the 
relativistic effects. 

Generally, it is believed that the ADAFs are hot and 
thick because of their inefficiency in radiating out the 
dissipated energy. In particular, as more energy is ad- 
verted, the disc becomes thicker. But our analysis may 
slightly change this picture in particular regarding to 
the thickness of the disc if the force due to the self- 
gravity of the disc in the vertical direction is considered. 
Our disc still rotates with nearly Keplerian profile, but 
its thickness is reduced significantly because of the self- 
gravity. In other words, we may have a fully advected 
disc, but its thickness is reduced in comparison to a 
case without self-gravity Also, as the disc becomes 
more advective, the effect of self-gravity becomes more 
evident. We think the present solutions are suitable for 
the outer parts of an ADAF not in inner parts. 

It is difficult the evaluated the precise picture of Ad- 
vection Dominated Accretion Flows (ADAFs) in the 
presence of B-field and self-gravity with self-similar 
method. However, this method can reproduce overall 
dynamical structure of the disc whit a set of given phys- 
ical parameter. Although our preliminary self-similar 
solutions are too simplified, they clearly improve our 
understanding of physics of ADAFs around a black 
hole. 
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